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Abstract. We explain how event-by-event fluctuations of particle ratios can 
constrain and falsify the statistical model of particle production in heavy ion 
collisions, using Kj-n fluctuations as an example. We define an observable 
capable of determining which statistical model, if any, governs freeze-out in 
ultrarelativistic heavy ion collisions. We calculate this observable for Kj-n 
fluctuations, and show that it should be the same for RHIG and LHC energies, 
as well as independent of centrality, if the Grand-Ganonical statistical model is 
an appropriate description and chemical equilibrium applies. We describe what 
happens in case of deviations from this scenario, such as light quark chemical 
non-equilibrium, strange quark over-saturation and local conservation (canonical 
ensemble) for strange quarks. We also introduce a similar observable capable, 
together with the published K* /K measurement, of ascertaining if an interacting 
hadron gas phase governs the system between thermal and chemical freeze-out, 
and of ascertaining its duration and impact on hadronic chemistry 
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One of the main objectives of heavy ion physics is to study the collective properties 
of strongly interacting matter. It's equation of state, degree of equilibration and 
phase structure, and the dependence of these on energy and system size. The natural 
approach to study particle production in heavy ion collisions is therefore statistical 
mechanics [U [H [3l |4] . In the last years, a consensus has developed that the statistical 
hadronization model can indeed fit most or all particles yields for AGS,SPS and RHIC 
energies [3 [6l H E [9] . 

This consensus, is, however, somewhat superficial. While it is true that one can 
get a reasonably nice-looking fit with a statistical model ansatz, it does not follow that 
the system is actually thermally and chemically equilibrated at freeze-out: Considering 
the paucity of data points when particle abundances are modeled, such a fit is by 
itself not a guarantee of the physical significance of parameters such as temperature 
and chemical potential. When statistical significance of these fits is calculated, it 
is apparent that the statistical model is nowhere near "proven" according to the 
standards generally accepted in particle physics. Particle fluctuations are a promising 
observable to falsify the statistical model and to constrain its parameters (choice of 
ensemble, strangeness/light quark chemical equilibrium) [lOj . One can immediately 
see that fluctuations are a stringent statistical model test by considering the fluctuation 
of a ratio between two random variables. 
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Since, for an equilibrated system, ((AA^i)^) ~ (N) ~ (V), where (V) is the system 
volume [llj . it is clear that cr'^^^j^^ depends on the hadronization volume in a manner 
opposite to that of particle yields, inversely rather than directly linearly proportionally. 
Volume fluctuations (which make a comparison of statistical model calculations to 
experimental data problematic), both resulting from dynamics and from fluctuations 
in collision geometry, should not alter this very basic result since volume cancels out 
event by event |12| . provided hadronization volume is the same for all particles (a 
basic statistical model requirement). 

Thus, observables such as (A^i) provided (A^i.2) and o'^^/jv2 measured 

using the same kinematic cuts, should be strictly independent of multiplicity and 
centrality, as long as the statistical model holds and the physically appropriate 
ensemble is Grand Canonical. If the temperature and chemical potentials between 
two energy regimes are approximately the same at freezeout (this should be the case 
for RHIC top energies and LHC, provided chemical equilibrium holds), this observable 
should also be identical. This could be used as a stringent test of the statistical model. 

Fluctuations are more sensitive to acceptance cuts than yields. A partial "fix" 
for acceptance cuts that does not require detector-specific analysis is mixed event 
subtraction, based on the idea that fluctuation effects resulting from acceptance cuts 
are present both in real and mixed events (this is the case for fluctuations, but not for 
correlations). Thus, an appropriate observable to model would be 

'^dyn — — fornix (2) 

where cr^„jj, is the mixed event width. In the absence of any correlations, this reduces to 
the Poisson expectation. (T^y„ (Known as [T3j '^fj"/pf^ when measured using histograms 
rather than event-by-event) is currently subject of intense experimental investigation 
[HI [T3] . We therefore propose to use 

= (3) 
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to test the statistical model among different energy, system size and centralities. 

"Primordial" fluctuations of each observable, ((ATVi ^)^), are calculable from 
Textbook methods [llj . but must also include corrections from resonance decays. 
Similarly, the correlation term in Eq. [l]is also given, in the Grand Canonical ensemble, 
by resonance abundance 

(A7ViAiV2)-^B,^i2(7V,) (4) 
j 

SHAREv2.X [161 can calculate all ingredients of any ^j^J/jy, ' incorporating the 
effect of all resonance decays, as well as chemical (non) equilibrium. 

It is important to underline that the value of 'I'j^J/jv2 (calculated from statistical 
model parameters and, as we will see, sensitive to the degree of chemical equilibration 
of the system) should be constant across any system where the intensive parameters 
are the same, for example different centrality regimes or system sizes at the same 
energy. For instance, since the chemical potential of Cu-Cu 200 GeV collisions should 
be comparable the chemical potential at Au-Au, 



A /TT 



dyn 

Au—Au ' 



{^7r)\cu-Cu dyn 



(5) 

Au-Au 



Cu-Cu (A^.)U„_^„ 

A quantitative test for equilibrium statistical hadronizationcan also be made 
between RHIC and LHC energies. Equilibrium thermal and chemical parameters 
are very similar at RHIC and the LHC (the baryo-chemical potential will be lower at 
the LHC, but, since both baryochemical potentials are small, the difference should 
not significantly affect tt and K abundance). Thus, ^'^-y^- should be identical, to 
within experimental error, for both the LHC and RHIC, over all multiplicities were 
the statistical model is thought to apply. 

According to |TS], chemical conditions at freeze-out (at SPS, RHIC and LHC) 
deviate from equilibrium, and reflect the higher entropy content and strangeness 
per entropy of the early deconfined phase through an over-saturated phase space 
occupancy for the light and strange quarks (7s > 7, > 1). If this is true, than 
^nI/N2 should still be independent of centrality for a given energy range, but should 
go markedly up for the LHC from RHIC, because of the increase in 7^ and 7s. 

We have calculated ^j^J/jy^ for RHIC and LHC energies, for the sets of parameters 
used in [18]. The left and right panel in Fig. [1] shows what effect three different sets 
of 7g,s inferred in [18] would have on "^^^-/j^- and In the left panel we 

have also included the value of ^'^--y^^- for top energy RHIC. As shown in [T^], this 
value for top centrality matches expectations for non-equilibrium freeze-out (and is 
significantly above equilibrium freeze-out). A centrality scan of ^'J--/^- , necessary to 
confirm the consistency of this result has not, however, as yet been performed. 

If non-statistical processes (mini-jets, string breaking etc.) dominate event-by- 
event physics, the flat ^j^]^/^^ scaling on centrality/multiplicity should be broken, and 

^Ari/Af2 '^oul'i exhibit a non-trivial dependence on Npart or dN/dy. 

This is also true if global correlations persist, such as is the case if the Canonical 
and micro-canonical ensembles [50] are physically more appropriate to describe the 
system. Hence, if strangeness at RHIC/the LHC is created and maintained locally, 
^iVi/Af2 should develop a "wiggle" at low centrality, and be considerably lower than 
Grand Canonical expectation. For 'i'^+zx- should be lower by a factor of two. 
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Figure 1. (color online) ^'T,, 



(Left panel) and 



(right panel) 



calculated for various statistical hadronization parameters |18| at the LHC. The 
left panel also shows the RHIC calculation |19l 



Fluctuations can also be useful in ascertaining the degree of reinteraction after 
an (assumed) statistical hadronization. A still unresolved ambiguity of statistical 
models is the duration, and impact on hadronic observables, of the phase between 
hadronization (the moment at which particles become the effective degrees of freedom) 
and freeze-out (the moment at which particles stop interacting). 

If, as commonly thought, chemical freezeout temperature is Tchem ^170 MeV 
and thermal freezeout temperature is Ttherm ~ 100 MeV, it follows that there is 
a significant interacting hadron gas phase that has the potential of altering all soft 
hadronic signatures. The failure to solve problems associated with HBT [3T], however, 
suggests that something fundamental is still not understood about freeze-out, and 
more direct probes of freeze-out dynamics are needed. 

The measurement of Resonance yields offers such a probe |22j . since short-lived 
hadronic resonances decay before the interacting hadron gas phase (if it exists) is over. 
Thus, rescattering of decay products can deplete the amount of observable resonances, 
while regeneration could create additional resonances not present at hadronization. 

The observation A(1520) and K*{%92) [23], at abundances below equilibrium 
statistical model expectations, could be interpreted as an indication of such 
reinteraction, with rescattering predictably dominating over regeneration. This 
interpretation, however, is not unique: Chemical non-equilibrium fits recover the 
resonance abundance exactly, with no need for an interacting hadron gas phase [22! ■ 

Thus detection of short lived resonances can not tell us freezeout dynamics unless 
a different signature, more sensitive to chemical freeze-out, is obtained. As is apparent 
[24l[T2] from Eqs[T]and|4]the correlation term is precisely the required observable, since 
correlations between multiplicities are fixed at chemical freeze-out. This correlation 
term can be measured by comparing observables such as *^+/.^- (correlated by 

-ftr*(892)) with (not correlated by resonances). We therefore define 

(i^*(892) ^ K+TT- 
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this result is somewhat spoiled by finite baryochemical potential, as well as higher 
lying resonances. To ascertain the size of these corrections, we have used SHARE 
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Figure 2. (color online) ANl/^y^ and K* /K calculated within the equilibrium 
and non-equilibrium statistical models 



to calculate both the K* /K and A'l'J.y^. As shown in Fig. [2] (left panel), these 
corrections make up a 10 % effect, less than the expected experimental error and not 
enough to alter the difference between a single freeze-out and two simultaneous ones. 

A long reinteracting hadron gas phase would in general bring the observed 
(final) abundance of K* / K away from the chemical freezeout value (either up, by 
regeneration, or down, by rescattering) . Thus, A^E*^^^ would become different from 
K* /K^ (Fig- m right panel). In the weak interaction limit, regeneration would 
presumably be rarer than rescattering so the observable K* / K~ abundance would 
be suppressed by a factor that combines the interaction width AF = ^rescattering — 
^regeneration with the lifetime of the interacting phase r 

1^ ^Avp-;^exp[-AFr] (7) 

observed 

In the strong reinteraction limit, rescattering and regeneration would re-equilibrate 
to a lower freeze-out Ttherm, so the observed K* / would be sensitive to the 
temperature difference as well as the mass difference (Am) between K* and K 



K* 



-Avf-^/^exp 

observed 



Am 



(8) 



One untested effect that could spoil this result is strong rescattering capable of bringing 
the K and tt out of the detector's acceptance region in phase space [inilll]. This effect 
can not be taken into account by mixed event techniques described in the previous 
section, and calculating it in a model-independent way is problematic. 

Inferring the presence of such a correction is however relatively straight-forward 
(Fig. [llright panel): The probability of such rescattering strongly depends on the 
width of the acceptance region. Thus, if A^PJ-^^ as a function of the rapidity window 
should go from zero (at small rapidity no multiplicity correlation survives) and saturate 
at a constant value (where the full resonance derived correlation is recovered). This 
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constant value, as long the rapidity window is much smaller than the total extent of 
rapidity of the system, is the quantity that can be related to 

chem 

The dependence on centrality of A^*^,^^, on the other hand, has to remain flat, 
since in the Grand Canonical limit the ratio of two particles should be independent of 
centrality, and the total system size should not alter the probability of a local process 
(such as scattering in/out of the acceptance region) to occur. 

If A^'J.^^ obeys the consistency checks described here it should be taken as an 

indication that the A\l/^^^ measurement reflects the value of at chemical freezeout. 

If A^J.^^ depends on rapidity up to the acceptance region of the detector, a more 
thorough effort to account for acceptance corrections to the correlation is needed. This 
can be done by using the same techniques utilised for direct resonance reconstruction 
[25] . Such an endeavour is detector specific, and goes beyond this write-up. 

In conclusion, we have shown that observables incorporating both yields and 
fluctuations give a stringent test of statistical models. We have also argued that such 
observables can be used to infer the duration of the interacting hadron gas phase, 
and its effect on hadronic observables. We expect that forthcoming experimental 
data, together with the methods elucidated here, will allow us to clarify some of the 
outstanding puzzles apparent in the study of heavy ion collisions. 
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